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I. Introduction 

This report describes an empirical analysis of a motion 
model that has been used to generate random submarine tracks for an 
ASW tactical decision aid. The model describes a submarine's 
motion as a series of transitions between the square cells of a 
grid that covers a defined operating region. A 3 X 3 transition 
matrix is associated with each cell of the grid which determines 
the submarine's transitions from a cell. The set of transition 
matrices define a Markov process. Despite its discrete nature, 
this Markov track generating process has been called a diffusion 
process in antisubmarine warfare tactical decision aid literature. 
The transition matrices are determined by tracks generated by an 
auxiliary stochastic process that is presumed to be of higher 
fidelity but more costly to implement than the Markov process. 

The auxiliary track generating process that was used in the 
analysis is a random tour process that was developed to generate 
operationally realizable submarine tracks. The process, which is 
described in detail in Reference 1, generates random tracks that 
have specified end points and a specified length or, for a 
submarine with a constant track speed, a specified transit time. 
The auxiliary process is implemented by the program that is listed 
in Appendix 1. 

The auxiliary tracks are the basis for the generation of two 
types of diffusion transition matrices by a procedure described in 
Section II. The first type, called static transition matrices, are 
the transition matrices that define the subject Markov motion 
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model. A static transition matrix is associated with each cell of 



the operating region grid. The second type, called dynamic 
transition matrices, define an alternate Markov motion model. 

This model differs from the subject model in that the transition 
matrix associated with a cell depends on the number of transitions 
that have taken place in the motion, that is, it is time dependent. 

The auxiliary tracks are also the basis for the generation of 
a set of maps that describe a submarine's position at a sequence of 
equally spaced times that are determined by a fixed time step. The 
number of times a cell is occupied at the end of a time step 
divided by the number of tracks is an estimate of the probability 
that the submarine was within the cell at the time step. For this 
reason, such maps are called probability maps in the following 
discussion . 

Tracks generated by the subject static diffusion Markov 
process and those generated by the alternate dynamic diffusion 
Markov process also can be used to generate probability maps. Some 
differences between the three track generating processes is evident 
by a comparison of the random tour probability maps, the dynamic 
diffusion probability maps and the static diffusion probability 
maps that are shown on the following pages. 
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II. The Transition Matrix Generation Procedure 



The transition matrices that are used in the analysis were 
determined by the program listed in Appendix 1. In the program, a 
random tour track is generated by a submarine moving with constant 
speed between two points located in a plane. The points lie on the 
x-axis of a rectangular coordinate system whose origin is midway 
between the points. Lines parallel to the x-axis and lines 
parallel to the y-axis define a grid of square cells in the plane. 
(For examples of the program geometry, see Section V.) In the 
following discussion of the transition matrix determination 
procedure, grid cells are identified by the rectangular coordinates 
of their centers. 

During its motion, the cell occupied by the submarine is 
determined at a sequence of times separated by a fixed time step. 
From this determination, transition matrices are constructed as 
follows: First, the cell that is occupied is determined for each 
time step. For I = 1 to NS (the number of time steps) , if cell 
(LO, MO) is occupied by the submarine at the end of the I-l st time 
step, at the end of I th time step, the submarine will be in a cell 
(LO + N, MO + K) where N = -1, 0, or 1 and K = -1, 0, or 1 . In 
this case, 1 is added to the (N,K) element of a 3 X 3 transition 
matrix associated with cell (L0,M0) for time 1-1. The elements of 
this matrix are stored in an array MT(N, K, LO, MO, 1-1) . In 
addition, 1 is added to the (N,K) element of a second 3X3 
transition matrix associated with cell (L0,M0). And, the elements 
of this matrix are stored in an array MTS(N, K, LO, MO). By 
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repeating this procedure for each random tour track, and dividing 
the resulting matrix elements by the number of tracks, the elements 
of the 3X3 matrix in the array MT become estimates of the 
elements of the dynamic transition matrix, and the elements of the 
matrix of the 3X3 matrix in the array MTS become estimates of 
the elements of the static transition matrix. 

The transition matrices and the probability maps used in the 
analysis were generated from 10,000 random tour tracks. Since the 
time step duration is equal to the cell side length divided by the 
track generating submarine's speed; during a time step, the 
submarine will either transition to an adjacent cell or remain in 
its current cell. This results in a transition matrix of nine 
elements, each element giving the probability of a transition. 
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III. Some Comparisons of the Track Generating Processes 

The dynamic and static transition matrices of the two discrete 
Markov track generating processes generate tracks that are random 
but differ in fine structure from those generated by the random 
tour process. This is not surprising considering the loss of 
information that occurs in their generation. Since the Markov 
motion is from a cell center to the cell center or a neighboring 
cell center, only nine motions are possible. In particular, for a 
submarine that moves with constant track speed, this results in a 
loss of position and velocity information. Relative to the 
auxiliary track generating process, the diffusion tracks no longer 
satisfy the constraints of constant speed, specified length and 
specified final end point. 

The fidelity of the diffusion tracks could be improved in two 
ways: First, by reducing the duration of the time step and 
consequently the size of the cells, less positional information 
would be lost. Second, by generating a velocity (direction of 
motion) distribution for each cell and each time step, transition 
matrices could be constructed that would be more descriptive of the 
random tour motion and less directional information would be lost. 
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IV. Conclusions 



The primary conclusion from the analysis is: For discrete 
Markov processes whose state space if current location and time, 
there are tactically reasonable submarine motion scenarios that can 
not be closely approximated by tracks generated with static 
transition matrices but can be closely approximated with dynamic 
transition matrices. There are, however, other tactically 
reasonable submarine motion scenarios that can not be closely 
approximated with either procedure. In what is possibly the 
simplest example of this, consider a submarine that moves among 
three adjacent cells labeled 1, 2 and 3. Suppose the submarine 
chooses either path (1,2,3) or path (3,2,1), each with probability 
.5. Given the submarine is in cell 2, the probability that the 
submarine will transition to cell 3 is dependent on its previous 
location as well as on its current location. Consequently, a 
Markov process whose state space is only current position and time 
can not describe its motion. 

For the random tour scenarios that were investigated, dynamic 
diffusion tracks were found to rather closely emulate the random 
tour tracks from which the dynamic transition matrices were 
derived. In addition, and in contrast to the static diffusion 
tracks, the dynamic diffusion tracks generated probability maps 
that duplicated the random tour probability maps. Although the 
static tracks were generally similar in appearance to the dynamic 
diffusion tracks, they did not uniformly terminate at the random 
tour track final end point. However, for the chosen scenario, 
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approximately half of the tracks did do so as is evident from the 
static diffusion probability maps that are shown in Section V. 

The divergence between the dynamic diffusion probability maps 
and the static diffusion probability maps is even more evident when 
the initial point and final point of a track coincide. The 
explanation for this is the variation of the dynamic transition 
matrices with time. That is, the same cell must perform a 
"submarine expansion" function early in the motion and a "submarine 
contraction" function towards the end. But, when the starting and 
ending cells are distant from each other, the function of any cell 
does not change significantly with time. That is, cells near the 
start cell are for always for "expansion" and those near the end 
cell are always for "contraction." 

The random tour process on which the empirical analysis is 
based was chosen to amplify the divergence between the random tour 
tracks, the dynamic diffusion tracks and the static diffusion 
tracks. The significance of this divergence for a tactical 
decision aid incorporating a motion model based on static 
transition matrices will depend on the desired applications of the 
decision aid. For example, it might be less significant if the 
location of the final end point of a submarine track were chosen 
randomly or if computer memory or processing speed limited modeling 
choices . 
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The figures in Section V that follows show probability maps 
and tracks generated by the random tour process, the dynamic Markov 
process transition matrices, and the static Markov process 
transition matrices. These figures illustrate the relative 
effectiveness of the dynamic and static Markov motion model in 
emulating the random tour motion model. The diffusion tracks that 
are shown were generated with 17 diffusions. This results in 
tracks that are approximately equal in length to the length the 
random tour tracks. 



8 



V. The Analysis Program Output 

The figures, tables and maps that follow are based on data 
generated by the analysis program that is listed in Appendix 1. 
They provide the basis for the empirical analysis. 
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Figure 1 . A track generated by a random tour process used in the 
analysis. The larger circles on the left and the right are the 
initial and final point of the track. The smaller circles are 
track positions at step times. 
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Figure 2. A superposition of 10,000 random tour tracks. These 
tracks determined the Markovian transition matrices that 
generated the diffusion probability maps that follow. 
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RANTRACK.BAS random tour probability map 
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random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 
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RANTRACK . BAS random tour probability map 
number of time steps = 4 



0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


2 


2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


3 


16 


42 


54 


22 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


8 


45 


176 


272 


309 


20 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


3 


22 


126 


337 


722 


832 


294 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


6 


31 


177 


449 


949 


1175 


485 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


5 


29 


126 


377 


746 


859 


271 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


17 


61 


137 


307 


315 


29 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


15 


29 


61 


32 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 
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RANTRACK.BAS random tour probability map 



number of time steps = 5 
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submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cel l entry sum = 10000 
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RANTRACK.BAS random tour probability map 
number of time steps = 6 
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submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 735 1 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 
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RANTRACK . BAS random tour probability map 
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steps = 
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submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cel l entry sum = 10000 
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RANTRACK . BAS random tour probability map 



nunber of time steps = 8 
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submarine speed in knots = 10 

distance between end points in nautical 
length of a cell side in nautical miles 
time step in hours = .3 
maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles - 5 

leg length distribution index = 3 

random number generator = standard 
random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 
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RANTRACK.BAS random tour probability map 



number of time steps = 9 
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submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 
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submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random nunber seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 
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submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 
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leg length distribution index = 3 
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number of tracks completed = 10000 

number of tracks terminated = 0 

cel l entry sum = 10000 
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submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 
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submarine speed in knots - 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 
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RANTRACK . BAS random tour probability map 
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RANTRACK.BAS random tour probability map 



nunber of 
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distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cel l entry sum = 10000 
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Table I. A dynamic diffusion transition matrix that was used to 
generate the dynamic diffusion map for 1 transition. 



RANTRACK.BAS dynamic transition matrix 

transition number = 1 

x = -5 
y = 0 



87 


1161 


961 


316 


1561 


3733 


80 


1133 


968 



submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = . 3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 
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Figure 3 A track generated by a dynamic process. The dynamic 
transition matrices were determined by the random tour tracks 
shown in Figure 2. The number of diffusions is 17. The circles 
are defined in Figure 1. 
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Figure 4 A superposition of 10,000 dynamic diffusion tracks. 
Their dynamic transition matrices were determined by the random 
tour tracks shown in Figure 2. The number of diffusions is 17. 
The circles are defined in Figure 1. 
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RANTRACK . BAS dynamic diffusion probability map 
number of diffusions = 0 
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RANTRACK.BAS dynamic diffusion probability map 
number of diffusions = 1 



0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 



0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

87 1161 961 0 

316 1561 3733 0 

80 1133 968 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



submarine speed in knots = 10 

distance between end points in nautical miles = 

length of a cell side in nautical miles = 3 

time step in hours = .3 
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delta in nautical miles = 5 

leg length distribution index = 3 
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number of tracks completed = 10000 

number of tracks terminated = 0 
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RANTRACK.BAS dynamic diffusion probability map 
number of diffusions = 2 
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distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 



35 



RANTRACK. BAS dynamic diffusion probability map 
number of diffusions = 3 
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submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 
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RANTRACK . BAS dynamic diffusion probability map 
number of diffusions = 4 
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RANTRACK.BAS dynamic diffusion probability map 
number of diffusions = 5 
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RANTRACK. BAS dynamic diffusion probability map 
number of diffusions = 6 
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RANTRACK.BAS dynamic diffusion probability map 



ntinber of diffusions = 7 
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RANTRACK.BAS dynamic diffusion probability map 



number of diffusions = 8 
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RANTRACK.BAS dynamic diffusion probability map 



number of diffusions = 9 
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RANT RACK - BAS dynamic diffusion probability map 



number of diffusions = 10 
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RANTRACK . BAS dynamic diffusion probability map 



number of diffusions = 11 
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RANTRACK . BAS dynamic diffusion probability map 



number of diffusions = 12 
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RANTRACK . BAS dynamic diffusion probability map 
number of diffusions = 13 
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number of tracks terminated = 0 

cell entry sum = 10000 
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RANTRACK . BAS dynamic diffusion probability map 
number of diffusions = 14 
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leg length distribution index = 3 

random nur vber generator = standard 
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number of tracks completed = 10000 

number of tracks terminated = 0 
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RANTRACK . BAS dynamic diffusion probability map 
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RANTRACK.BAS dynamic diffusion probability map 
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RANTRACK.BAS dynamic diffusion probability map 
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submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 
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Table II. A static diffusion transition matrix that was used to 
generate the static diffusion maps. 



RANTRACK . BAS 


static 


transition matrix 




X = -5 
y = 0 


107 


1407 


1290 




375 


2008 


4848 




99 


1389 


1332 



submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 



51 




Figure 5 A track generated by a static diffusion process. The 
static transition matrices were determined by the random tour 
tracks shown in Figure 2. The number of diffusions is 17. The 
circles are defined in Figure 1. 
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Figure 6 A superposition of 10,000 static diffusion tracks. 
Their static transition matrices were determined by the random 
tour tracks shown in Figure 2. The number of diffusions is 17. 
The circles are defined in Figure 1. 
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RANTRACK. BAS static diffusion probability map 
nurber of diffusions = 0 
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submarine speed in knots = 10 

distance between end points in nautical miles 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maxi mun nurber of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random n tuber generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 

cell entries are rounded integer 



RANTRACK.BAS static diffusion probability map 
number of diffusions = 1 
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distance between end points in nautical miles = 30 

length of a cell siae in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random nuntoer generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

nun±>er of tracks terminated = 0 

cell entry sum = 10000 

cell entries are rounded integer 
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RANTRACK.BAS static diffusion probability map 
number of diffusions = 2 
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submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximun rumber of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 

cell entries are rounded integer 
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RANTRACK.BAS static diffusion probaoility map 
number of diffusions = 3 
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submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 10000 

cell entries are rounded integer 
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RANTRACK.BAS static diffusion probability map 
number of diffusions = 4 
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submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random n inter generator = standard 

random n under seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sun = 10000 

cell entries are rounded integer 
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RANTRACK. BAS static diffusion probability map 
number of diffusions = 5 
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RANTRACK. BAS static diffusion probability map 
nunber of diffusions = 6 
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RANTRACK.BAS static diffusion probability map 



number of diffusions = 7 
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RANTRACK. BAS static diffusion probability map 



number of diffusions = 8 
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RANTRACK.BAS static diffusion probability map 



number of diffusions = 9 
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RANTRACK . BAS static diffusion probability map 



number of diffusions = 10 



0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


1 


2 


5 


0 


0 


0 


0 


1 


3 


10 


22 


0 


0 


0 


1 


3 


9 


26 


58 


0 


0 


0 


1 


5 


16 


44 


110 


0 


0 


0 


1 


5 


17 


50 


122 


0 


0 


0 


2 


6 


17 


46 


109 


0 


0 


0 


1 


3 


9 


26 


58 


0 


0 


0 


0 


1 


3 


10 


25 


0 


0 


0 


0 


0 


1 


3 


8 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



submarine speed in knots = 10 

distance between end points in nautical 

length of a cell side in nautical miles 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles » 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sum = 9999.999 

cell entries are rounded integer 



0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


8 


5 


2 


0 


0 


0 


0 


0 


0 


52 


34 


14 


3 


0 


0 


0 


0 


0 


205 


179 


93 


21 


2 


0 


0 


0 


0 


461 


461 


319 


131 


19 


0 


0 


0 


0 


612 


701 


673 


448 


136 


0 


0 


0 


0 


453 


452 


315 


125 


18 


0 


0 


0 


0 


206 


178 


87 


20 


2 


0 


0 


0 


0 


57 


35 


14 


3 


0 


0 


0 


0 


0 


8 


5 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



30 



0 0 

0 0 

0 0 

0 0 

0 0 

1 0 

8 10 

40 51 

117 174 

210 354 

252 420 

213 341 

116 176 

41 56 

7 9 

1 1 

0 0 

0 0 

0 0 

0 0 

0 0 

mi les = 

= 3 



64 



RANTRACK.BAS static diffusion probability map 



nunber of diffusions = 11 
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RANTRACK. BAS static diffusion probability map 



number of diffusions = 12 
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RANTRACK.BAS static diffusion probability map 



number of diffusions = 13 
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length of a cell side in nautical miles = 3 
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RANTRACK.BAS static diffusion probability map 



nunber of diffusions = 14 
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RANTRACK.BAS static diffusion probability map 



number 
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RANTRACK.BAS static diffusion probability map 



number of diffusions = 16 
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RANT RACK . BAS static diffusion probability map 



number of diffusions = 17 
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24 


41 


52 
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10 


24 


50 


95 


162 


263 


606 
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0 


0 


2 


5 


12 


28 


64 


131 


277 


930 
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0 


0 


0 



submarine speed in knots = 10 

distance between end points in nautical miles = 30 

length of a cell side in nautical miles = 3 

time step in hours = .3 

maximum number of time steps = 17 

track length in nautical miles = 50 

delta in nautical miles = 5 

leg length distribution index = 3 

random number generator = standard 

random number seed = 7351 

number of tracks completed = 10000 

number of tracks terminated = 0 

cell entry sun = 10000 

cell entries are rounded integer 
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VI. The Analysis Program 

The analysis program requires the following inputs: the speed 
of the submarine in knots; the distance between the end points in 
nautical miles; the track length in nautical miles; delta, the 
average distance between course changes (leg length) in nautical 
miles; the leg length distribution parameter, an integer that 
defines the gamma distribution that determines leg length (Inputing 
1 gives an exponential distribution.) and the number of tracks to 
be generated. If an auxiliary pseudorandom number generator is not 
chosen, the option to supply a seed for the QuickBASIC random 
number generator is given. Next, the option to either view single 
tracks and terminate the program or to input the number of tracks 
to be generated to produce data of the kind presented above. 

The program is written for a computer monitor with screen mode 
12 capability. To change this requirement, the statement SCREEN 12 
on Line 4780 must be changed. For example, for a monitor with CGA 
capability, it could be replaced by SCREEN 2. 

To run the program under QuickBASIC, the /ah switch must be 
used to duplicate the program results in this report. This switch 
allows arrays of size greater than 64K. In addition, it may be 
necessary to have FILES = 8 and BUFFERS =2 in the CONFIG.SYS 
file and to remove TSR file activation statements from the 
AUTOEXEC.BAT file. 

The program code listed in Appendix 1 provides for the 
addition of an auxiliary random number generator as an alternative 
to the QuickBASIC generator through Lines 200, 210 and 220. These 



72 



lines allow auxiliary generator setup code in a subroutine starting 
at Line 7720 and generator code in a subroutine starting at Line 
7810. This line number was chosen to accommodate code that is 
based on a Generalized Feedback Shift Register (GFSR) pseudorandom 
number generator that is listed in Reference 3. If an auxiliary 
generator will not be added. Lines 200, 210 and 220 can be removed. 
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Appendix 1. The Analysis Program Listing 



10 NS = "RANTRACK.BAS": REM a program to analyze some random track generation procedures 
20 CLS : DEFINT I-N 

30 PI = 4 * ATN(1 ) : NTR = 0: RS = "standard": FLAG3 = 0: FLAG4 = 0: FLAGS = 0 
40 PRINT : INPUT "submarine speed in knots"; V 

50 PRINT : INPUT "distance between the end points in nautical miles"; D 
60 PRINT : INPUT "track length in nautical miles"; TL 
70 IF TL < D THEN 60 

80 DM = SQR(TL * TL • D * 0) / 2: REM maximum distance from the x-axis 

90 SL = DM / 7: IF D / 2 > SL THEN SL = (D / 2) / INT (D / 2 / SL): REM length of a cell side 

100 TS ■ SL / V: REM time step in hours 
110 MNS = 30: REM maximum number of steps 

120 NS = CINTCTL / SL): REM neuter of time steps in the track 
130 IF NS > MNS THEN GOTO 60 

140 SF = (64 / 35) * DM / 1000: REM screen factor 

150 AS = "cell dimension = " + STRS(SL) ♦ " nm": B$ = "time step =" ♦ STRS(TS) + » hr" 

160 PRINT : PRINT AS: PRINT : PRINT BS 

170 PRINT : INPUT "delta in nautical miles"; DEL: REM maximun average leg length in nautical miles 

180 IF DEL < TL * TL / (TL - D) / 50 THEN 170 

190 PRINT : INPUT "leg length distribution index"; NI 

200 AS = " ": PRINT : INPUT "standard or auxiliary random number generator (s/a)"; AS 
210 IF AS = "A" OR AS = "a" THEN GOSUB 7720: GOTO 280 

220 IF AS * "S" OR AS = "s" THEN 230 ELSE 200 

230 AS = " ": PRINT : INPUT "supply a random number seed (y/n)"; AS 

240 IF AS * "N" OR AS = "n" THEN 280 

250 IF AS = "Y" OR AS = "y» THEN FLAG4 = 1 ELSE 230 

260 PRINT : INPUT "random number seed"; RNS 

270 IF RNS < -32768! OR RNS > 32767 THEN 260 ELSE RANDOMIZE RNS 

280 AS = " ": PRINT : INPUT "view only single tracks (y/n)"; AS 

290 IF AS = "Y" OR AS = "y« THEN FLAG5 = 1: NT = 1: GOTO 320 

300 IF AS = "N» OR AS = "n" THEN 310 ELSE 280 

310 PRINT : INPUT "number of tracks"; NT 

320 A0 = TL / 2: REM semimajor axis of the initial focusing ellipse 

330 CO = D / 2: REM distance from the center to a focus of the initial focusing ellipse 
340 MNL = 30 

350 LL = 9: LU = 9: ML = 10: MU = 10: MT = 6 
360 REDIM MP(-LL TO LU, -ML TO MU, MNS) 

370 REDIM MT ( - 1 TO 1, -1 TO 1, -LL TO LU, -ML TO MU, MNS) 

380 REDIM MTS( - 1 TO 1, -1 TO 1, -LL TO LU, -ML TO MU) 

390 FOR J = 1 TO NT 

400 REDIM X(MNL + 1), Y(MNL + 1), XS(MNL + 1), YS(MNL + 1), FR(MNL + 1), T (MNL + 1), TLEG(MNL + 1) 

410 X(0) = -D / 2: Y(0) = 0: REM coordinates of the initial point for all tracks 

420 FR(0) = 0: TLEG(0) = 0 

430 FLAG1 = 0 

440 A = A0: C = CO: S = 0 

450 GOSUB 4770 

460 FOR I = 1 TO MNL: REM leg number 
470 B1 = SQR( (TL ♦ D) / (TL - D)) 

480 B2 = SQR( 1 - D * D / TL / TL) 

490 GOSUB 7700 

500 ON ERROR GOTO 510: GOTO 520 
510 RESUME 490 
520 FLAG2 = 0 

530 IF RAND <= .5 THEN 550 

540 RAND = RAND - .5: FLAG2 = 1 

550 RC = 2 * ATN(B2 * TAN(PI * RAND - ATN(BD) + D / TL): REM course in radians relative to the minor axis of 

the focusing ellipse 

560 ON ERROR GOTO 0 

570 IF FLAG2 = 1 THEN RC = PI - RC 

580 AL = (TL - D) / (TL - D * SIN(RC)): : REM average leg length for the relative course RC 

590 GOSUB 7620 

600 ON ERROR GOTO 0 

610 FR ( I ) = RC ♦ S: REM course in radians on leg I 

620 TLEGM = (A * A - C * C) / (A - C * SIN(RC) ) : REM maximum leg length in nautical miles for the relative 

course RC 

630 IF TLEG >= TLEGM THEN TLEG = TLEGM: FLAGl = 1: REM track is terminated at the final point 
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640 GOSUB 6030 

650 LINE (X(I), Y( I ) )* (X( I - 1), Y(I - 1)) 

660 A = Air REM remaining track length in nautical miles 

670 C = Cl: REM distance between the final point and the end point of the previous leg 

680 S = BR - PI / 2: REM direction of the minor axis of the focusing ellipse 

690 TLEG(I) * TLEG: REM leg length of leg I 

700 IF FLAG1 = 1 THEN 750 

710 NEXT I 

720 NTR = NTR + 1 

730 J = J - 1 

740 GOTO 1280: REM track terminated 
750 LINE (D / 2, 0>-(X(I), Y< I )) 

760 NLEG = I ♦ 1: REM number of legs 

770 TLEG = 2 * C: TLEG(NLEG) = TLEG 

780 H * S + PI / 2: REM bearing in radians of final point from the end point of the leg 

790 FR(NLEG) = H: REM course in radians on last leg 

800 X(NLEG) = D / 2: REM x-coordinate of final point for all tracks 

810 Y(NLEG) = 0: REM y-coordinate of final point for all tracks 

820 T ( 0 ) = 0 

830 FOR I = 1 TO NLEG 

840 T * TLEG(I) / V: REM time in hours to traverse leg I 

850 T(I) = T(I - 1) ♦ T: REM time in hours from the initial point to the end point of leg I 
860 NEXT I 

870 XS(0) = -D / 2: YS(0) = 0 
880 K = 1: TSK = TS 
890 FOR I = 1 TO NLEG 

900 IF TSX > TO) AND I < NLEG THEN 960: REM the next time step point is past the Ith leg 

910 TK = TSX - TO - 1): REM time from the starting point of the Ith leg to the Xth time step point 

920 XS(X) = X(I - 1) ♦ V * TK * SIN(FRO)): REM x-coordinate of the Xth time step point 

930 YS(K) = Y(I - 1) ♦ V * TX * COS(FRO)): REM y-coordinate of the Xth time step point 

940 X = X ♦ 1: TSX « TSX ♦ TS: REM next time step 

950 IF X > NS THEN 980 ELSE 900 

960 NEXT I 

970 ERASE X, Y, FR, T, TLEG 
980 FOR X = 1 TO NS - 1 
990 CIRCLE (XS(X), YS(K>), 8 * SF 
1000 NEXT X 

1010 IF FLAGS = 0 THEN 1050 

1020 AS INPUT “Quit (y/n)“; AS 

1030 IF AS = "N" OR AS = “n" THEN CLS : GOTO 400 

1040 IF AS = “Y" OR AS = »y» THEN END ELSE 1020 

1050 FOR I = 0 TO NS 
1060 FOR L = -LL TO LU 

1070 IF (L - 1 / 2) * SL < XS(I> AND XS(!) <= (L + 1 / 2) * SL THEN 1090: REM the x-coordinate of the Ith time 
step point is in range L 
1080 NEXT L 

1090 FOR M = -ML TO MU 

1100 IF (M - 1 / 2) * SL < YS( I ) AND YS(I) <= (M + 1 / 2) * SL THEN 1130: REM the y-coordinate of the Ith time 
step point is in range M 
1110 NEXT M 
1120 GOTO 1270 

1130 MP(L, M, I) = MP(L, M, I) + 1: REM random tour probability map element 

1140 IF I > 0 THEN 1170 

1150 L0 = L: M0 = M 

1160 GOTO 1270 

1170 LI = L: Ml = M 

1180 FOR N = -1 TO 1 

1190 FOR X = -1 TO 1 

1200 IF NOT (Li = L0 + N AND Ml = M0 + X) THEN 1240 

1210 MT(N, X, L0, M0, I - 1) = MT(N, X, L0, M0, I - 1) + 1: REM dynamic diffusion matrix element 

1220 MTS(N, X, L0, M0) = MTS(N, X, L0, M0) + 1: REM static diffusion matrix element 

1230 GOTO 1260 

1240 NEXT X 

1250 NEXT N 

1260 L0 = Li: M0 = Ml 

1270 NEXT I 

1280 NEXT J 

1290 AS = " “: PRINT : INPUT “display a random tour probability map (y/n)"; AS 
1300 IF AS = “N" OR AS = »n» THEN 1500 



75 



1310 IF AS = »Y» OR A$ = "y" THEN 1320 ELSE 1290 
1320 SCREEN 0: CLS 
1330 PRINT "NS = NS 

1340 PRINT : INPUT "enter the nunber of time steps from 0 to NS"; I 
1350 IF I < 0 OR I > NS THEN 1340 
1360 FOR M = 10 TO -10 STEP -1 

1370 XS = STR$(MP( -9, M, I)): X$ = LTR IM$(RTR IMS(XS) ) : PRINT XS; TAB (5); 

1380 FOR N = -8 TO 8 

1390 XS = STRS(MP(N, M, I)): XS = LTR1M$(RTRIMS(XS)): PRINT XS; TAB((N ♦ 9) * 4 ♦ 5); 
1400 NEXT N 

1410 XS = STRS(MP(9, M, I)): XS = ITRIHS(RTRIM$(XS)): PRINT XS 
1420 NEXT M 



1430 AS = » ": PRINT : INPUT 
1440 IF AS = "N" OR AS = "n" 

1450 IF AS = "V" OR AS = "y" 

1460 WIDTH IPRINT 130 
1470 GOSUB 4880 
1480 WIDTH IPRINT 80 
1490 GOTO 1290 

1500 AS = " PRINT : INPUT 
1510 IF AS = »N" OR AS = »n" 

1520 IF AS = »Y" OR AS = "y" 



"print the random tour probability map (y/n)"; AS 
THEN 1290 

THEN 1460 ELSE 1430 



"print a random tour probability map (y/n)“; AS 
THEN 1610 

THEN 1530 ELSE 1500 



1530 SCREEN 0: CLS 
1540 PRINT "NS = "; NS 

1550 PRINT : INPUT "enter the nunber of time steps from 0 to NS I 

1560 IF I < 0 OR I > NS THEN 1550 

1570 WIDTH LPRINT 130 

1580 GOSUB 4880 

1590 WIDTH LPRINT 80 

1600 GOTO 1290 

1610 AS = " ": PRINT : INPUT "print all the random tour probability maps (y/n)"; AS 

1620 IF AS = "N" OR AS = "n" THEN 1680 

1630 IF AS * "Y" OR AS = "y" THEN 1640 ELSE 1610 

1640 WIDTH LPRINT 130 

1650 FOR 1 = 0 TO NS 

1660 GOSUB 4880 

1670 NEXT I 



1680 AS = " ": PRINT : INPUT "print all the random tour probability maps to a file (y/n)" 

1690 IF AS = "N" OR AS = "n" THEN 1800 

1700 IF AS = "Y" OR AS = "y" THEN 1710 ELSE 1680 

1710 PRINT : INPUT "input the data file name"; FS 

1720 ON ERROR GOTO 1730: GOTO 1740 

1730 RESUME 1710 

1740 OPEN "O", #1, FS 

1750 FOR I = 0 TO NS 

1760 GOSUB 7050 

1770 NEXT I 

1780 CLOSE #1 

1790 ON ERROR GOTO 0 

1800 ERASE HP 



1810 AS = " ": PRINT : INPUT "display a dynamic transition matrix (y/n)"; AS 

1820 IF AS = "N" OR AS = "n" THEN 1920 

1830 IF AS = "Y" OR AS = "y" THEN 1840 ELSE 1800 

1840 SCREEN 0: CLS : GOSUB 6480 

1850 PRINT : PRINT 

1860 FOR K = 1 TO -1 STEP -1 

1870 PRINT : PRINT MT(-1, K, L, M, 1 - 1); TAB( 15); MT(0, K, L, M, I - 1); TAB(29); MT(1, 
1880 NEXT K 

1890 PRINT : PRINT : PRINT : PRINT "nunber of tracks completed = NT 
1900 PRINT : PRINT "nunber of tracks terminated = "; NTR 
1910 PRINT : GOTO 1810 

1920 AS = " PRINT : INPUT "print a dynamic transition matrix (y/n)"; AS 

1930 IF AS = "N" OR AS = "n" THEN 2100 

1940 IF AS = "Y" OR AS = "y" THEN 1950 ELSE 1920 

1950 SCREEN 0: CLS : GOSUB 6480 

1960 AS = NS + " dynamic transition matrix": LPRINT AS 
1970 LPRINT : LPRINT "transition nunber = "; I 
1980 LPRINT : LPRINT "x = "; L 
1990 LPRINT "y = "; M 



; as 



K. 1, 



M, 



1 ) 
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2000 FOR K = 1 TO -1 STEP -1 

2010 LPRINT : IPRINT ; TAB(15); MT(-1, X, L, H, I - 1); TA8(29); MT(0, X, L, M, I - 1); TAB(43); MT(1, X, L 
M, I - 1) 

2020 NEXT K 

2030 IPRINT : IPRINT 

2040 GOSUB 6130 

2050 IF FLAG4 = 0 THEN N = 34 ELSE N = 33 

2060 FOR K = 1 TO N 

2070 LPRINT 

2080 NEXT K 

2090 GOTO 1920 

2100 AS = " PRINT : INPUT "print nonzero dynamic transition matrices (y/n)"; AS 

2110 IF AS = "N" OR AS = "n" THEN 2390 

2120 IF AS = "Y" OR AS = "y" THEN 2130 ELSE 2100 

2130 SCREEN 0: CLS 

2140 PRINT "NS = "; NS 

2150 PRINT : INPUT "enter the transition number from 1 to NS"; I 
2160 IF I < 1 OR I > NS THEN 2150 

2170 AS = NS + " nonzero dynamic transition matrices": LPRINT AS 
2180 LPRINT : LPRINT "transition nutter = I: LPRINT 
2190 FOR L = -(LL - 1) TO (LU - 1) 

2200 FOR H = -(ML - 1) TO (MU - 1) 

2210 MAT = 0 

2220 FOR N = -1 TO 1 

2230 FOR X = -1 TO 1 

2240 MAT = HT(N, X, L, M, I - 1) + MAT 

2250 NEXT K 

2260 NEXT N 

2270 IF MAT * 0 THEN 2340 
2280 LPRINT : LPRINT "x = "; L 
2290 LPRINT "y = M 
2300 FOR K = 1 TO -1 STEP -1 

2310 LPRINT : LPRINT ; TA8(15); MT ( - 1 , X, L, M, I - 1); TAB(29); MT(0, X, L, M, I - 1); TAB (43); MT(1, X, L 
M, I - 1) 

2320 NEXT X 

2330 LPRINT : LPRINT 

2340 NEXT M 

2350 NEXT L 

2360 LPRINT : LPRINT 

2370 G0SU8 6130 

2380 GOTO 2100 

2390 AS = " PRINT : INPUT "print a dynamic transition matrix to a file (y/n)"; AS 

2400 IF AS = "N" OR AS = "n" THEN 2630 

2410 IF AS = "Y" OR AS = "y» THEN 2420 ELSE 2390 

2420 SCREEN 0: CLS : INPUT "input the data file name"; FS 

2430 ON ERROR GOTO 2440: GOTO 2450 

2440 RESUME 2420 

2450 OPEN "0", #1, FS 

2460 GOSUB 6480 

2470 AS = NS + " dynamic transition matrix": PRINT #1, AS 
2480 PRINT #1, : PRINT #1, "transition number * "; I 
2490 PRINT #1, : PRINT #1, "x = L 
2500 PRINT #1, "y = »; M 
2510 FOR X = 1 TO -1 STEP -1 

2520 PRINT #1, : PRINT #1, ; TAB( 15); MT(-1, X, L, M, I - 1); TAB(29); MT(0, X, L, M, I - 1); TAB(43); MT(1 
X, L, M, I - 1) 

2530 NEXT X 

2540 PRINT #1, : PRINT #1, 

2550 GOSU8 6560 

2560 IF FLAG4 = 0 THEN N = 27 ELSE N = 26 
2570 FOR X = 1 TO N 
2580 PRINT #1, 

2590 NEXT X 
2600 CLOSE #1 
2610 ON ERROR GOTO 0 
2620 GOTO 2390 

2630 AS = " ": PRINT : INPUT "generate dynamic diffusion tracks (y/n)"; AS 

2640 IF AS = "Y» OR AS = "y“ THEN 2660 

2650 IF AS = "N" OR AS = "n" THEN 2770 ELSE 2630 
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2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 



SCREEN 0: CLS 



PRINT "NS = ", NS 

PRINT : INPUT "enter the number of diffusions from 1 to NS"; ND 

IF ND < 1 OR ND > NS THEN 2680 

PRINT : INPUT "enter the number of tracks"; NOT 

CLS 

GOSUB 4770 
FOR J = 1 TO NOT 
GOSUB 5050 
NEXT J 
GOTO 2630 

AS = *' PRINT : INPUT "generate a dynamic diffusion probability map (y/n)"; AS 

IF AS = "Y" OR AS = »y" THEN 2800 

IF AS = »N" OR AS = "n" THEN 3050 ELSE 2770 

SCREEN 0: CLS 



PRINT "NS = "; NS 

PRINT : INPUT "enter the nunber of diffusions from 0 to NS "; ND 

IF ND < 0 OR ND > NS THEN 2820 

REDIM MPD(-LL TO LU, -ML TO MU) AS INTEGER 

GOSUB 5330 

AS = " ": PRINT : INPUT "display the dynamic diffusion probability map (y/n)"; AS 

IF AS = "N» OR AS = "n" THEN 2970 

IF AS = "Y" OR AS = »y" THEN 2890 ELSE 2860 

SCREEN 0: CLS 

FOR M = 10 TO -10 STEP -1 

XS = STRS(MPD( -9, M)): XS = LTRIMS(RTRIMS(XS) ) : PRINT XS; TAB(5 ); 

FOR N = -8 TO 8 

XS = STRS(MPD (N , M)): XS = LTRIMS(RTRIMS(XS) ) : PRINT XS; TAB( (N ♦ 9) * 4 ♦ 5); 
NEXT N 

XS = STRS(MPD(9, M)): XS = LTRIMS(RTRIMS(XS) ) : PRINT XS 
NEXT M 



AS = » PRINT : 
IF AS = »N" OR AS 
IF AS = "Y" OR AS 
WIDTH LPRINT 130 
GOSUB 5680 
WIDTH LPRINT 80 
ERASE MPD 
GOTO 2770 
AS = " ": PRINT : 
IF AS = "N» OR AS 
IF AS = "Y» OR AS 



INPUT "print the dynamic diffusion probability map (y/n)"; AS 
= "n" THEN 2770 
= "y" THEN 3000 ELSE 2970 



INPUT "print all the dynamic diffusion probability maps (y/n)"; 
= "n" THEN 3160 
= "y" THEN 3080 ELSE 3050 



AS 



WIDTH LPRINT 130 
FOR ND = 0 TO NS 

REDIM MPD( -LL TO LU, -ML TO MU) AS INTEGER 

GOSUB 5330 

GOSUB 5680 

NEXT ND 

ERASE MPD 

WIDTH LPRINT 80 

AS = " ": PRINT : INPUT "print all the dynamic diffusion probability maps to a file (y/n)"; 

IF AS = "N" OR AS = "n" THEN 3310 

IF AS = "Y" OR AS = "y" THEN 3190 ELSE 3160 

PRINT : INPUT "input the data file name"; FS 

ON ERROR GOTO 3220 

GOTO 3230 

RESUME 3190 

OPEN "O", #1, FS 

FOR ND = 0 TO NS 

REDIM MPD( -LL TO LU, -ML TO MU) AS INTEGER 

GOSUB 5330 

GOSUB 6700 

NEXT ND 

CLOSE #1 

ON ERROR GOTO 0 

ERASE MPD 

AS = " ": PRINT : INPUT "display a static transition matrix (y/n)"; AS 

IF AS = "N" OR AS = "n" THEN 3430 

IF AS = "Y" OR AS = "y" THEN 3350 ELSE 3320 



AS 



78 



3350 SCREEN 0: CLS : GOSUB 6510 

3360 PRINT : PRINT 

3370 FOR X = 1 TO -1 STEP -1 

3380 PRINT : PRINT HTSC-1, X, l, M); TAB(15); MTS(0, X, L, M); TAB(29); MTS(1, X, L, M) 

3390 NEXT K 

3400 PRINT : PRINT : PRINT : PRINT "nurber of tracks completed = NT 

3410 PRINT : PRINT "nurber of tracks terminated = "; NTR 

3420 PRINT : GOTO 3320 

3430 AS = " PRINT : INPUT "print a static transition matrix (y/n)"; AS 

3440 IF AS = "N" OR AS = "n» THEN 3600 

3450 IF AS = "Y" OR AS = "y" THEN 3460 ELSE 3430 

3460 SCREEN 0: CLS : GOSUB 6510 

3470 AS = NS + “ static transition matrix": LPRINT AS 

3480 LPRINT : LPRINT "x = L 

3490 LPRINT "y = »; M 

3500 FOR X = 1 TO -1 STEP -1 

3510 LPRINT : LPRINT ; TAB(15); MTSC-1, K, L, M); TAB(29); MTS(0, X, L, M); TAB(43); MTS(1, X, L, M) 

3520 NEXT X 

3530 LPRINT : LPRINT 

3540 GOSUB 6130 

3550 IF FLAG4 = 0 THEN N = 36 ELSE N = 35 

3560 FOR X = 1 TO N 

3570 LPRINT 

3580 NEXT X 

3590 GOTO 3430 

3600 AS = " PRINT : INPUT "print all nonzero static transition matrices (y/n)"; AS 

3610 IF AS = "N" OR AS = "n" THEN 3840 

3620 IF AS = "Y" OR AS = "y" THEN 3630 ELSE 3600 

3630 SCREEN 0: CLS 

3640 AS = NS ♦ " nonzero static transition matrices": LPRINT AS 
3650 FOR L = -<LL - 1) TO (LU - 1) 

3660 FOR M = -(ML - 1) TO (MU - 1) 

3670 MAT = 0 

3680 FOR N = -1 TO 1 

3690 FOR X = -1 TO 1 

3700 MAT = MTS(N, X, L, M) ♦ MAT 

3710 NEXT X 

3720 NEXT N 

3730 IF MAT = 0 THEN 3800 
3740 LPRINT : LPRINT »x = "; L 
3750 LPRINT "y = "; M 
3760 FOR X = 1 TO -1 STEP -1 

3770 LPRINT : LPRINT ; TAB(15); MTS(-1, K, L, M); TAB(29); MTS(0, K, L, M); TAB(43); MTS(1, K, L, M) 

3780 NEXT K 

3790 LPRINT : LPRINT 

3800 NEXT M 

3810 NEXT L 

3820 LPRINT : LPRINT 

3830 GOSUB 6130 

3840 AS = " ": PRINT : INPUT "print a static transition matrix to a file (y/n)"; AS 

3850 IF AS = "N" OR AS = "n" THEN 4080 

3860 IF AS = »Y» OR AS = "y" THEN 3870 ELSE 3840 

3870 SCREEN 0: CLS : INPUT "input the data file name"; F$ 

3880 ON ERROR GOTO 3890: GOTO 3900 
3890 RESUME 3870 
3900 OPEN "O", #1, FS 
3910 GOSUB 6510 

3920 AS = NS + " static transition matrix" 

3930 PRINT #1, AS 

3940 PRINT #1, : PRINT #1, "x = "; L 
3950 PRINT #1, "y = "; M 
3960 FOR X = 1 TO -1 STEP -1 

3970 PRINT #1, : PRINT #1, ; TAB(15); MTS(-1, K, L, M); TAB(29); MTS(0, X, L, M); TA8(43); MTS(1, X, L, M) 
3980 NEXT X 

3990 PRINT #1, : PRINT #1, 

4000 GOSUB 6560 

4010 IF FLAG4 = 0 THEN N = 30 ELSE N = 29 
4020 FOR X = 1 TO N 
4030 PRINT #1, 
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4040 

4050 

4060 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 

4720 



NEXT K 
CLOSE #1 
ON ERROR GOTO 0 
GOTO 3850 

AS = " PRINT : INPUT "generate static diffusion tracks (y/n)"; A$ 

IF AS = "Y" OR AS = "y" THEN 4110 

IF AS = "N" OR AS = "n» THEN 4230 ELSE 4080 

SCREEN 0: CLS 



PRINT "NS = ", NS 

PRINT : INPUT "enter the rumber of diffusions from 1 to NS"; ND 

IF ND < 1 OR ND > NS THEN 4130 

PRINT : INPUT "enter the number of tracks"; NOT 

CLS 

GOSUB 4770 
FLAG3 * 1 
FOR J = 1 TO NOT 
GOSUB 5050 
NEXT J „ 

GOTO 4080 

AS = " ": PRINT : INPUT "generate a static diffusion probability map (y/n)"; AS 

IF AS = "Y" OR AS = "y" THEN 4260 

IF AS = "N" OR AS = "n" THEN 4520 ELSE 4230 

SCREEN 0: CLS 

PRINT "NS = "; NS 

PRINT : INPUT "enter the number of diffusions from 0 to NS "; ND 
IF ND < 0 OR ND > NS THEN 4280 
RED IH PDT(-LL TO LU, -ML TO MU) 

GOSUB 7220 
GOSUB 7550 

AS = " ": PRINT : INPUT "display the static diffusion probability map (y/n)"; AS 
IF AS = "N" OR AS = "n" THEN 4440 



IF AS = "Y" OR AS = "y" THEN 4360 ELSE 4330 

SCREEN 0: CLS 

FOR M = 10 TO -10 STEP -1 

XS = STRS(CI NT (PDT ( -9, M))): XS = LTRIMS(RTRIMS(XS)) : PRINT XS; TAB(5); 

FOR N = -8 TO 8 

XS = STR$(CINT(PDT(N, M))): XS = LTR IMS(RTRIMS(X$) ) : PRINT XS; TAB( (N + 9) * 4 + 5); 
NEXT N 

XS = STRS(CINT (PDT (9, M))): XS = LTRIMS(RTRIMS(XS) ) : PRINT XS 
NEXT M 



AS = " PRINT : 
IF AS = "N" OR AS 
IF AS = "Y" OR AS 
WIDTH LPRINT 130 
GOSUB 5850 
WIDTH LPRINT 80 
ERASE PDT 
GOTO 4230 
AS = « ": PRINT : 
IF AS = "N" OR AS 
IF AS = "Y" OR AS 



INPUT "print the static diffusion probability map (y/n)"; AS 
= "n" THEN 4230 
= "y" THEN 4470 ELSE 4440 



INPUT "print all the static diffusion probability maps (y/n)"; AS 
= "n" THEN 4630 
= "y» THEN 4550 ELSE 4520 



WIDTH LPRINT 130 
FOR ND = 0 TO NS 
REDIM PDT ( -LL TO LU, -ML TO MU) 

GOSUB 7220 
GOSUB 5850 
NEXT ND 
ERASE PDT 
WIDTH LPRINT 80 

AS = " ": PRINT : INPUT "print all the static diffusion probability maps to a file (y/n)"; 

IF AS = "N» OR AS = "n» THEN 4760 

IF AS = "Y" OR AS = "y" THEN 4660 ELSE 4630 

PRINT : INPUT "input the data file name"; FS 

ON ERROR GOTO 4680: GOTO 4690 

RESUME 4660 

OPEN "O", #1, FS 

FOR ND = 0 TO NS 

REDIM PDT( -LL TO LU, -ML TO MU) 

GOSUB 7220 



AS 
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4730 GOSUB 6870 
4740 NEXT ND 
4750 ERASE PDT 
4760 END 

4770 XM = 10 * SL: YM = 8 * SL 
4780 SCREEN 12 

4790 WINDOW ( -XM, -YM)-(XM, YM) 

4800 CIRCLE (-D / 2, 0), 10 * SF: CIRCLE (D / 2, 0), 10 * SF 
4810 FOR M = 0 TO 19 

4820 LINE ( -XM + (M ♦ 1 / 2) * SL, YM - SL / 2)-<-XM ♦ (M ♦ 1 / 2) * SL, -YM + SL / 2) 

4830 NEXT M 

4840 FOR M = 0 TO 15 

4850 LINE (-XM + SL / 2, YM - (M + 1 / 2) * SL)-(XM - SL / 2, YM - (M + 1 / 2) * SL) 

4860 NEXT M 
4870 RETURN 

4880 At = Nt + " random tour probability map": LPRINT At 
4890 LPRINT : LPRINT "number of time steps = I: LPRINT 
4900 FOR M = 10 TO -10 STEP -1 

4910 Xt = STRt(MP( -9, M, I)): Xt = LTRIM!(RTRIMS(XS)): LPRINT Xt; TAB(MT); 

4920 FOR N = -8 TO 8 

4930 Xt = STRt(MP(N, M, I)): Xt = LTRIMt(RTRIMt(Xt)): LPRINT Xt; TAB((N ♦ 9) * (MT - 1) + MT); 
4940 NEXT N 

4950 Xt = STRt(MP(9, M, I)): LPRINT Xt 
4960 NEXT M 
4970 GOSUB 6130 
4980 GOSUB 6270 

4990 LPRINT "cell entry sum = "; SM 

5000 IF FLAG4 = 0 THEN N = 22 ELSE N = 21 

5010 FOR X = 1 TO N 

5020 LPRINT 

5030 NEXT K 

5040 RETURN 

5050 LO = -D / 2 / SL: MO = 0: REM initial cell indices 
5060 FOR I = 1 TO ND 

5070 REDIM MTAC-1 TO 1 , -1 TO 1) AS INTEGER 

5080 MPSUM = 0 

5090 FOR N = -1 TO 1 

5100 FOR X = -1 TO 1 

5110 IF FLAG3 = 0 THEN 5130 

5120 MPSUM = MTS(N, X, LO, MO) + MPSUM: GOTO 5140 

5130 MPSUM = MT(N, K, LO, MO, I - 1 ) ♦ MPSUM 

5140 MTA(N , K) = MPSUM 

5150 NEXT X 

5160 NEXT N 

5170 GOSUB 7700 

5180 RNU = RAND * MPSUM 

5190 FOR N = -1 TO 1 

5200 FOR X = -1 TO 1 

5210 IF RNU < MTALN, X) THEN 5260: REM transition from (LO+N.MO+X) to (LO.MO) 

5220 NEXT X 
5230 NEXT N 
5240 ERASE MTA 
5250 GOTO 5320 

5260 LI = LO + N: Ml = MO + X 

5270 LINE (LI * SL, Ml * SL)-(LO * SL, MO * SL) 

5280 CIRCLE (LI * SL, Ml * SL), 8 * SF 

5290 LO = LI: MO = Ml 

5300 ERASE MTA 

5310 NEXT I 

5320 RETURN 

5330 REDIM MSUM( -LL TO LU, -ML TO MU, MNS) 

5340 FOR I = 0 TO NS - 1 
5350 FOR L = -LL TO LU 
5360 FOR M = -ML TO MU 
5370 FOR N = -1 TO 1 
5380 FOR X = -1 TO 1 

5390 MSUM(L , M, I) = MT (N, X, L, M, I) + MSUM(L, M, I) 

5400 NEXT X 
5410 NEXT N 
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5420 NEXT M 
5430 NEXT L 
5440 NEXT I 

5450 MPD(-D / 2 / SL, 0) = NT: REM initial cell probability factor 

5460 IF NO = 0 THEN 5660 

5470 FOR I = 1 TO NO 

5480 REDIM B(-ll TO LU, -ML TO MU) 

5490 FOR L = -(LL - 1) TO (LU - 1) 

5500 FOR M ® -(ML - 1) TO (MU - 1) 

5510 FOR N = -1 TO 1 
5520 FOR K = -1 TO 1 

5530 IF MSUM( L - N, M - X, I - 1 ) = 0 THEN 5550 

5540 B(L, M) = MT(N, X, L - N, M - X, I - 1 ) / MSUM(L - N, M - X, I - 1) * MPD(L - N, M - X) + B(L, M) 

5550 NEXT X 
5560 NEXT N 
5570 NEXT M 
5580 NEXT L 

5590 FOR L * -LL TO LU 

5600 FOR M = -ML TO MU 

5610 MPD(L, M) = B(L, M): REM dynamic diffusion probability map element 

5620 NEXT M 

5630 NEXT L 

5640 ERASE B 

5650 NEXT I 

5660 ERASE MSUM 

5670 RETURN 

5680 AS = NS + " dynamic diffusion probability map": LPRINT A$ 

5690 LPRINT : LPRINT "nunber of diffusions = ND: LPRINT 
5700 FOR M = 10 TO -10 STEP -1 

5710 XS = STR$(HPD( -9, M)): X$ = LTRIM$(RTRIM$(X$) ) : LPRINT XS; TAB(MT); 

5720 FOR N = -8 TO 8 

5730 XS * STR$(MPD(N, M)): X$ = LTRIM$(RTRIM$(X$) ) : LPRINT XS; TAB( (N + 9) * (MT - 1) + MT); 

5740 NEXT N 

5750 XS = STR$(MPD(9, M)): XS = LTRIM$(RTRIM$(X$) ) : LPRINT XS 
5760 NEXT M 
5770 GOSUB 6130 
5780 GOSUB 6340 

5790 LPRINT "cell entry sum = SM 

5800 IF FLAG4 = 0 THEN N = 22 ELSE N = 21 

5810 FOR X = 1 TO N 

5820 LPRINT 

5830 NEXT X 

5840 RETURN 

5850 AS = NS + " static diffusion probability map": LPRINT AS 
5860 LPRINT : LPRINT "number of diffusions = ND: LPRINT 
5870 FOR M = 10 TO -10 STEP -1 

5880 XS = STRS(CINT(PDT(-9, M))): XS = LTRIM$(RTR1M$(X$)): LPRINT XS; TAB(MT); 

5890 FOR N = -8 TO 8 

5900 XS = STR$(CINT(PDT(N, M))): XS = LTRIMS(RTR1MS(XS)): LPRINT XS; TAB( (N + 9) * (MT - 1) + MT); 

5910 NEXT N 

5920 XS = STRS(CINT(PDT(9, M))): XS = LTRIMS(RTRIMS(XS)): LPRINT XS 
5930 NEXT M 
5940 GOSUB 6130 
5950 GOSUB 6410 

5960 LPRINT "cell entry sun = "; SM 

5970 LPRINT "cell entries are rounded integer" 

5980 IF FLAG4 = 0 THEN N = 21 ELSE N = 20 

5990 FOR X = 1 TO N 

6000 LPRINT 

6010 NEXT X 

6020 RETURN 

6030 X(I) = X(I - 1) + TLEG * SIN(FR(I)): REM x-coordinate of the end of leg I 

6040 Y(I) = Y(I - 1) ♦ TLEG * COS(FR(I)): REM y-coordinate of the end of leg I 

6050 X = D / 2 - X(I): Y = 0 - Y(I): REM determine the bearing and range of the final point from the end of leg 
I 

6060 RF = SQR(X * X + Y * Y): REM range of the final point from the end of leg I 

6070 Cl = RF / 2: REM distance of the center of the trial focusing ellipse from the final point 

6080 A1 = A - TLEG / 2: REM one-half the remaining track length 

6090 IF ABS(X / RF) = 1 THEN BB = PI / 2 * SGN(X) ELSE BB = ATN(X / RF / S0R(1 - X * X / RF / RF)) 
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6100 IF ABS( Y / RF) = 1 THEN BR = PI / 2 * (1 - SGN(Y)) ELSE BR = PI / 2 - ATN( Y / RF / SQR ( 1 - Y * Y 
RF)) 

6110 IF BB < 0 THEN BR = 2 * PI - BR: REM bearing in radians of the final point from the end point of 
6120 RETURN 



6130 

6140 

6150 

6160 

6170 

6180 

6190 

6200 

6210 

6220 

6230 

6240 

6250 

6260 

6270 

6280 

6290 

6300 

6310 

6320 

6330 

6340 

6350 

6360 

6370 

6380 

6390 

6400 

6410 

6420 

6430 

6440 

6450 

6460 

6470 

6480 

6490 

6500 

6510 

6520 

6530 

6540 

6550 

6560 

6570 

6580 

6590 

6600 

6610 

6620 

6630 

6640 

6650 

6660 

6670 

6680 

6690 

6700 

6710 

6720 

6730 

6740 

6750 

6760 

6770 



LPRINT : LPRINT ‘'submarine speed in knots = “; V 

LPRINT “distance between end points in nautical miles = C 

LPRINT “length of a cell side in nautical miles = SL 

LPRINT “time step in hours =“; TS 

LPRINT “maximum number of time steps = NS 

LPRINT “track length in nautical miles = TL 

LPRINT “delta in nautical miles = “; DEL 

LPRINT “leg length distribution index = NI 

LPRINT “random number generator = “; RS 

IF FLAG4 * 0 THEN 6240 

LPRINT “random number seed * “; RNS 

LPRINT “number of tracks completed = NT 

LPRINT “number of tracks terminated = NTR 

RETURN 

SM = 0 

FOR L = -LL TO LU 

FOR M = -ML TO MU 

SM = MP(L, M, I) ♦ SM 

NEXT M 

NEXT L 

RETURN 

SM = 0 

FOR L = -LL TO LU 

FOR M = -ML TO MU 

SM = MPD( L , M) ♦ SM 

NEXT M 

NEXT L 

RETURN 

SM = 0 

FOR L = -LL TO LU 
FOR M = -ML TO MU 
SM = PDT ( L , M) + SM 
NEXT M 
NEXT L 
RETURN 

PRINT “NS = “; NS 

PRINT : INPUT “enter the transition number from 1 to NS"; I 
1 OR I > NS THEN 6490 
INPUT “x-coordinate"; L 
-(LL - 1) OR L > (LU - 1) THEN 6510 
INPUT “y-coordinate"; M 
-(ML - 1) OR M > (MU - 1) THEN 6530 



IF I < 
PRINT 
IF L < 
PRINT 
IF M < 
RETURN 
PRINT #1, 
PRINT #1, 
PRINT #1, 
PRINT #1, 



TL 



: PRINT #1, “submarine speed in knots = “; V 
"distance between end points in nautical miles = 

"length of a cell side in nautical miles = SL 
“time step in hours =“; TS 
PRINT #1, “maximum number of time steps = “; NS 
PRINT #1, "track length in nautical miles = 1 
PRINT #1, "delta in nautical miles = “; DEL 
PRINT #1, “leg length distribution index = “; NI 
PRINT U 1, "random number generator = RS 
IF FLAG4 = 0 THEN 6670 
PRINT #1, “random number seed = “; RNS 
PRINT #1, “number of tracks completed = NT 
PRINT #1, “number of tracks terminated = NTR 
RETURN 

AS = NS + " dynamic diffusion probability map": 

PRINT #1, : PRINT #1, “nurber of diffusions = 

FOR M = 10 TO -10 STEP -1 

XS = STRS(MPD( -9, M)): X$ = LTRIMS(RTRIMS(XS) ) : 

FOR N = -8 TO 8 

XS = STR$(MPD( N , M)): X$ = LTRIMS(RTRIMS(XS) ) : PRINT #1, X$; TAB( (N ♦ 9) 
NEXT N 

X$ = STR$(MPD(9, M)): XS = LTRIM$(RTRIM$(X$) ) : PRINT #1, XS 



PRINT #1, AS 
; ND: PRINT #1, 

PRINT #1, XS; TAB(MT); 



(MT - 1) + MT); 



/ RF / 
leg I 
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6780 NEXT H 
6790 GOSUB 6560 
6800 GOSUB 6340 

6810 PRINT #1, "cell entry sum = "; SH 
6820 IF FLAG4 = 0 THEN N = 16 ELSE N = 15 
6830 FOR K = 1 TO N 
6840 PRINT #1, 

6850 NEXT K 
6860 RETURN 

6870 AS = NS + " static diffusion probability map": PRINT #1, AS 
6880 PRINT DM, : PRINT #1, "nunber of diffusions = ND: PRINT #1, 

6890 FOR H = 10 TO -10 STEP -1 

6900 XS = STRS(C!NT(PDT(-9, M))): XS = LTRIM$(RTR!MS(XS)): PRINT #1, XS; TAB(MT); 

6910 FOR N = -8 TO 8 

6920 XS = STRS(CI NT (P0T(N , M))): XS * LTRIM$(RTR!HS(XS)) : PRINT #1, XS; TAB((N + 9) * (MT - 1) + MT) 
6930 NEXT N 

6940 XS = STRSCCI NT ( PDT ( 9 , M))): XS = ITRIMS(RTR!MS(XS)): PRINT #1, XS 
6950 NEXT H 
6960 GOSUB 6560 
6970 GOSUB 6410 

6980 PRINT #1, "cell entry sun = SM 

6990 PRINT #1, "cell entries are rounded integer" 

7000 IF FLAG4 = 0 THEN N = 16 ELSE N = 15 
7010 FOR K = 1 TO N 
7020 PRINT #1, 

7030 NEXT X 
7040 RETURN 

7050 AS = N$ + " random tour probability map": PRINT #1, AS 
7060 PRINT #1, : PRINT #1, "number of time steps = "; I: PRINT #1, 

7070 FOR H = 10 TO -10 STEP -1 

7080 XS = STRS(MP( -9, H, I)): XS = LTRIM$(RTRIM$(X$>): PRINT #1, XS; TAB(MT); 

7090 FOR N = -8 TO 8 

7100 XS = STRS(MP(N, M, I)): XS = LTRIM$(RTRIM$(XS)): PRINT #1, XS; TAB((N + 9) * (MT - 1) + MT); 
7110 NEXT N 

7120 XS = STR$(MP(9, M, I)): XS = LTRIMS(RTRIMS(XS)): PRINT #1, XS 
7130 NEXT M 
7140 GOSUB 6560 
7150 GOSUB 6270 

7160 PRINT #1, "cell entry sum = »; SM 
7170 IF FLAG4 = 0 THEN N = 16 ELSE N = 15 
7180 FOR X = 1 TO N 
7190 PRINT #1, 

7200 NEXT X 
7210 RETURN 

7220 REDIM MSUMS(-LL TO LU, -ML TO MU) 

7230 FOR L = -LL TO LU 

7240 FOR M = -ML TO MU 

7250 FOR N = -1 TO 1 

7260 FOR X = -1 TO 1 

7270 MSUMSCL, M) = MTS(N, X, L, M) + MSUMSCL, M) 

7280 NEXT X 
7290 NEXT N 
7300 NEXT M 
7310 NEXT L 

7320 PDT(-D / 2 / SL, 0) = NT: REM initial cell probability factor 

7330 IF NO = 0 THEN RETURN 

7340 FOR I = 1 TO ND 

7350 REDIM B( -LL TO LU, -ML TO MU) 

7360 FOR L = -(LL - 1) TO (LU - 1) 

7370 FOR M = -(ML - 1) TO (HU - 1) 

7380 FOR N = -1 TO 1 

7390 FOR X = -1 TO 1 

7400 IF HSUHS(L • N, M - K) ■ 0 THEN 7420 

7410 B(L, M) = MTS(N, K, L - N, M - K) / MSUMS(L - N, M - X) * PDT (L - N, M - X) + B(L, M) 

7420 NEXT X 
7430 NEXT N 
7440 NEXT M 
7450 NEXT L 

7460 FOR L = -LL TO LU 
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7470 FOR M * -ML TO MU 

7480 PDT(L, M) = B(L, M): REM static diffusion probability map element 

7490 NEXT M 

7500 NEXT L 

7510 ERASE B 

7520 NEXT I 

7530 ERASE MSUMS 

7540 RETURN 

7550 SM = 0 

7560 FOR L = -LL TO LU 

7570 FOR M = -ML TO MU 

7580 SM * PDT(L, M) + SM 

7590 NEXT M 

7600 NEXT L 

7610 RETURN 

7620 TLEG = 0 

7630 FOR L = 1 TO NI 

7640 GOSUB 7700 

7650 ON ERROR GOTO 7660: GOTO 7670 
7660 RESUME 7640 

7670 TLEG = -DEL * AL * LOGO - RAND) / NI ♦ TLEG: REM trial leg length generator for the relative course RC 
7680 NEXT L 
7690 RETURN 

7700 IF RS = "standard" THEN RAND = RND ELSE GOSUB 7810: REM auxiliary random number generator subroutine branch 
7710 RETURN 

7720 R$ * "auxiliary": REM the first line of an auxiliary random number generator setup subroutine 
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